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I. INTRODUCTION 



Intracules are two-particle density distributions obtained from the spinless second-order reduced density matrix [I] 

92 \ r r ^j= J ( r i> r 2' r 3 , . . . , t-jv) * {r[,r' 2 , r 3 , . . . , r N ) dr 3 ■ ■ ■ dr N , (1) 

where ^(ri, r 2 , r 3 , . . . , r^r) is the iV-particle position wave function. Intracules are usually normalized to the number 
of particle pairs N(N - l)/2. 

The seminal intracule is the Position intracule 

P{u)=jp,[ r r +U ; r r +u )drdn u , (2) 

which was introduced long ago by Coulson and Neilson [2] to study correlation effects in the helium atom. In pj, 
u = ri — r%, u = |u| = ri2 and is the angular part of it. P{u) gives the probability density for finding two 
particles separated by a distance u and has been widely studied [3H2"T]. 
The corresponding Momentum intracule 22,, [23J is 

M(v) = f p 2 ( ^ V ^ ' r + q )e lqv drdqdudri V7 (3) 

where t> = pi — u = = P12 and f2„ is the angular part of v. M(v) gives the probability density for finding two 
particles moving with a relative momentum v. 

Starting with the Wigner distribution [24] [25] , one can construct a family of intracules [HI [16] , which provide 
two-electron position and/or momentum information. Within this family, the patriarch is the Omega intracule 

n(u,v,u) = j^ j P*(r + u + q \ r r + iy iqv Ku-O uv )drdqdn u dn v , (4) 

where u = 9 UV is the dynamical angle between the vector u and v. fl(u,v,u>) can be interpreted as the joint quasi- 
probability density for u, v and uj. The quasi prefix emphasizes that Cl(u,v,u>) is not a rigorous probability density 
and, indeed, it may take negative values [23]. Based on the observation of Rassolov 26 that both relative position 
and relative momentum are important to describe the correlation between pairs of electrons, and because the Omega 
intracule contains information on both quantities, Cl(u, v, uS) has been extensively used in Intracule Functional Theory 

(IFT) [BII3127H33]. 

Appropriate integrations [TH] reduce the Omega intracule to lower-order intracules such as P(u), M(v) and the 
Angle intracule [27] 125] 

T(u>) = / / £l(u,v,uj)dudv } (5) 
Jo Jo 

which provides information on the angle uj between u and v. A similar reduction yields the Dot intracule |17l 130) 

D(x) = / / v ' - ' dzdu. (6) 
J J x uzsmuj 

The variable x = u ■ v = uv cos u> combines information on the relative position and momentum of the particles, and 
it is easy to show that it gives the rate of change of u 2 , i.e. 

HI- 2 - < 7 > 

In this way, x sheds light on the motion of the electrons. For example, x = implies that the electrons are moving 
in such a way that their separation is constant. This could arise, for example, if they were in a circular orbit around 
their centre of mass. 

Although D(x) is usually a non-negative function and has proven useful for understanding electronic behaviour |27j 
and for estimating electron correlation energies in atomic and molecular systems [29[ [30] , its connection to the Omega 
suggests that it is not a rigorous probability density. However, in the following Section, we show how to derive the 
exact probability distribution of x. 
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II. THE POSMOM INTRACULE 



Wc define the Posmom intracule X(x) to be the exact probability density for the variable x = u ■ v. It is the 
two-particle version of the Posmom density S(s) where s = r ■ p [34TI36] and, as we have argued that s describes 
particle trajectories, we now propose that x likewise characterizes pair trajectories. 

The quantum mechanical operator 



,fi(| +P . Vp ) (8) 



2 

is known to be an unbounded self-adjoint operator |371 38: and its two-particle equivalent is 

x = -2ih(~ + u- V u ^j , (9) 

where V is the gradient operator. Both s and x correspond to quantum mechanical observables. 

Following the same approach used in [34] to obtain S(s) from s, one can show that the Posmom intracule can be 
expressed as the Fourier transform 

X(x) = — / X{k)e ikx dk (10) 
of the two-particle hyperbolic autocorrelation function 

x(k) = J P 2 ( r + r ekhu ; 1 1 ^ ( ( ^") drdu - 

This expression can be simplified after defining the intracule density matrix 

/ a f (U-u/2 , U-u'/2\ , TT 

Pu{u , «)=JP2{y + J /2 ; U + u / /2 )dU, (12) 

where U = r% + V2 is the extracule vector, and yields 

X(k) = j p u (e +kh u , e- kh u) du. (13) 

In an entirely analogous way, the Dot intracule can be expressed as the Fourier transform 

1 f°° ^ 

D{x) = — / D{k)e lkx dk (14) 

of the /-Dot function (30] 



and the latter can be reduced to 



D(k) = J p u ((1 + kh)u , (1 - kh)u) du. (16) 



Comparing ( 13 1 with ( 16 1 and the Taylor expansion of the exponential function 

e ±** = l±jfe/i+— ft a + ... (17) 

reveals that the probability density D(x) derived from the Wigner distribution is a first-order approximation to the 
exact density X(x). Thus, the quasi-intracule is correct to 0(h) and becomes exact in the classical limit H 0. 
Remarkably, one can construct the exact density from the approximate density using the mapping 



D(ta,nh(kh)) 
cosh 3 (kh) 



4 



TABLE I. One- and two-particle hyperbolic autocorrelation functions. 



Density 



Intracule 



Dot 

Posmom 
Relation 



SW(fc) = J Pi ((1 + kh/2)r , (1 - kh/2)r) dr 
S(k) = J Pl (t 

~ = g w (2tanh(fcfi/2)) 
( ' cosh 3 (fcft/2) 

= Sw(k) + 0(h 2 ) 



D(k) = p u ((1 + kh)u , (1 - kh)uj 



<Iu 



X(k) 



pu I e u , e u 



m = g(taah(M)) 
cosh 3 (fc/i) 

= b{k) + o{h 2 ) 



Table I gives the one- and two-particle hyperbolic autocorrelation functions, the corresponding first-order (Wigner) 
approximations, and the relations between them. In Table [TJ the one-particle density matrix is given by 



Pi 



P2 



, r 2 , r 2 , 



dr 2 



(19) 



If the wave function is expanded in one-electron functions, <f> a (r), the reduced two-particle density matrix becomes 



P2 



,r 2 



^ Pabcd&a 



'ibcd 



where P a bcd is a two-particle density matrix clement. In this case, (11) is given by 



x ( k ) = ^2 Pabcd [abcd]^ , 



(20) 



(21) 



'ibcd 



where we have introduced the two-particle hyperbolic autocorrelation integral [abcd]^. For example, if the basis 
functions are s-type Gaussians, we obtain (in atomic units) 



\ssss 



x 



e -«|r-A| 2 e -^|r+smh(fc) U -B| 2 e -7|r+o X p(/c) M -C| 2 e -«5|r+co S li(fc)"-0| 2 dr du 
3 



J 3 / 2 



exp 



1 /li? 12 



where 



and 



C = a + /3 + 7 + <5, 

J = sinh 2 (k) + 7exp 2 (fc) + <5cosh 2 (fc)] - [/3sinh(fc) + 7exp(£;) + Scosh(k)] 2 , 
F = a[3\A - B\ 2 + aj\A - C\ 2 + aS\A - D\ 2 + (3j\B - C\ 2 + /3S\B - D\ 2 + 1 5\C - D\ 2 , 



H = smh(k)G B + exp(fc)G c + cosh(fc)G D , 
G B = a/3(B -A) + /3 7 (B - C) + (5S(B - D), 
G c = ai {C -A) + /3 7 (C - B) + 7 «5(C - D), 
G D = aS{D -A)+ (35(D - B) + jS(D - C). 



(22) 



(23a) 
(23b) 
(23c) 



(24a) 
(24b) 
(24c) 
(24d) 



The scalars and vectors above are independent of the choice of origin and X(k) and X(x) are therefore likewise 
independent. This contrasts with the annoying origin-dependence |35j of the one-electron posmom density S(s). 

Integrals of higher angular momentum can be generated by differentiating [ssss] % with respect to the Cartesian 
coordinates of the basis function centers, as first suggested by Boys [39], or, more efficiently, using recurrence relations 
[31j . We have written a program to compute X(x) within an spd Gaussian basis set and implemented this in a 
development version of the Q-Chem 3.2 quantum chemistry package [3D]. 
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Eqs (21 1 - p4| are easily modified to generate P(u), M(v), T(w) and D(k). In particular, if the functions sinh(fc), 
exp(fc) and cosh(fc) are replaced by their first-order approximations [k, 1 + k and 1, respectively) in the expressions 
for J and if, one obtains the /-Dot integrals [ssss]g [3D]. 

In the special case of concentric s-type Gaussians, the intracule integrals become 

r i 471-5/2 2 ( A* 2 \ , 
[ssssjp = , u exp I — —u , (25a) 



4^5/2 



JM-^372-«'exp^--^j, (25b) 

7T 3 (A — 2?7 cos 2 w) 

y sinw, (25c) 



2C 3 / 2 (A + 77 cos 2 w) 

^3 



£-3/2' 

^ 3 
J37 2 ' 



(25d) 
(25e) 



where 



M = (a + /?)( 7 + <5), (26a) 

i/ = (a + 7)C9 + <5), (26b) 

C = (a + 5)09 + 7), (26c) 

X = 4:{af3-f + a(3S + ajS + fij5), (26d) 



a + <5 ,5 + 7/ V" + 5 /? + 7 



^a + <5 ^ + 7 

if = e [Pk 2 + 7 (1 + fc) 2 + 6] - [/3k + 7(1 + k) + 8} 2 . (26g) 



In the calculations described below, we have computed X(x) and D(x) numerically using Eqs. (10) and (14 1. D(k), 
X(k), D(x) and X(x) are all even functions and we will therefore focus only on x > and k > 0. 

Physical interpretations of the variables u, v, ui and x are summarized in Fig. [T] The three limiting configurations 
lj = 0, 7r/2 and n (which correspond to x = u v, and —u v) are depicted for the weak (u and v large), medium (where 
one of u and v is large and the other is small) and strong correlation (u and v small) regimes. A faithful description 
of electron correlation requires information about the relative position u and momentum v, but also on the mutual 
orientation u) of these two vectors, which gives insight into the nature of the electrons' mutual orbit. The Dot and 
Posmom intracules provide information about the distribution of values of x — uvcosuj, and thus about the type of 
correlation regime (weak, medium or strong). However, as noted above, being a first-order approximation of X(x), 
the information gathered in D(x) is slightly biased. The effects of this approximation will be investigated below. 



In Section III the Posmom intracule is investigated alongside D(x), T(w), P(u) and M(v) for the two electrons in 
a parabolic quantum dot. In Section |IV[ we turn our attention to the electrons in a helium atom or helium-like ion. 
We also compare the Posmom intracules for ground and excited states and study the effect of the dimensionality of 
the space T>. Atomic units are used throughout. 



III. PARABOLIC QUANTUM DOTS 



In our study of the Posmom intracule in parabolic quantum dots [41j , we consider three different treatments of the 
Coulomb interaction between the two electrons. First, the non-interacting case, in which it is simply ignored; second, 
the Hartree-Fock (HF) case [12] in which it is approximated in a mean-field sense; third, the exact treatment which 
is possible for certain values of the harmonic confinement force constant [33] H3] . 



«_ 03=0 
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i co=n/2 
I j x =0 

co=n 

x =-u v (medium) 
medium correlation 


co=0 
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x =-u v (medium) 
medium correlation 
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FIG. 1. Physical interpretation of the variables u, v, uj and x in the weak, medium and strong correlation reg 



A. Hamiltonian and wave functions 

The Hamiltonian is 

n = - i(v? + v|) + v( ri ) + v(r 2 ) + — , 

2 r i2 

where 

2 

V M = h 

is the external harmonic potential and 1/k 2 is the force constant. 

The S ground state of the non-interacting system has the wave function 

*o(n,r 2 ) = tpo{ri)ipo(r 2 ), 
V^o(r) = (nn) 3 ^ 4 exp 



2 k 



and the energy 

3 



E = 

K 



The more accurate HF wave function 



*HF(n,r 2 ) = ^HF^l^HF^) 
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TABLE II. Energies of parabolic quantum dots for different treatments of the interelectronic interaction. 





N G 


K = 2 


K = 10 






1.5 


0.3 


Ehf 


1 

2 
3 
4 
5 

O'Neill and Giljj] 
RagolQ 


2.04 

2.038 439 
2.038 438 9 
2.038 438 8718 
2.038 438 871755 
2.038 438 87 
2.038 438 87176 


0.53 
0.529 04 
0.529 041 5 
0.529 041525 6 
0.529 041525 56 


£cxact 




2. 


0.5 



a Reference | 14| : 7 basis functions. 
b Reference |46| : 11 basis functions. 



is not known in closed form, but can be efficiently treated numerically by expanding i/'hfM in a Gaussian basis 

JVg 

^hf(t-) = ^2 c.j exp(-a i r 2 ). (33) 
i=i 

The HF energy can be directly minimized with respect to the coefficients Cj and exponents aj using a numerical solver 
[35] , thus avoiding the self-consistent field procedure usually needed for this kind of calculation [351 S3 ■ 
The exact wave function and energy can be found in closed form [44) for certain values of k. For example, for k = 2 

* 2 (ri, r 2 ) = (l + T -f) Vo(riM)(r 2 ), (34a) 

E 2 = 2, (34b) 

and, for n = 10, 

*io(n,r 2 ) =( 1 + r ~f + %) Mn)Mr2), (35a) 

E 10 = 1/2. (35b) 
Table [Tl] shows the convergence of Ehf with Nq for k = 2 and k = 10. The correlation energy 

.E c = -Ecxact — -Bhf (36) 
for k = 2 (J^c = 38.438 871 755 mE h ) agrees with earlier work [H US]- For k = 10, we find E c = 29.041 525 56 m£i, 



B. Position Intracule 



The non-interacting Position intracule is 



Po(u) 



- { u 2 exp 



u 

2k 



(37) 



and the HF intracule Phf,k(w) is found from (25a I and (33). For n = 2 and k = 10, the exact intracules are given by 



P2(U 

Pio(u) = 



8 + 5Jtt 



u exp 



u 



(1 + w/2 + w 2 /20) 2 , 
5/2(240 + 61v/5tt) " 



u exp 



u 
20 



(38) 
(39) 



Equation ( 38 ) has been reported previously [13] . 



C. Momentum Intracule 



Using the same notation as above, the non-interacting Momentum intracule is 

MoO) = ^Vexp ^-^J-) > 



Mhp,2(u) and A/hf, 10(f) are obtained from (25b) and (33), and the exact Momentum intracules are 

-1 2 



M 2 (v) 



8v' 



8 + 5^ 



1 

iv 



1 r I erf ( —= ) e 



V2j 



M w (v) 



48V5 + 6iy7r 



— + (4- 5v 2 ) e"T" 



where erf(z) is the error function [48 . Equation (41) has been reported previously [13] . 



(40) 

(41) 
(42) 



D. Angle Intracule 

The Angle intracule of two non-interacting particles is entirely determined by the Jacobian factor and is 

T oM = 2 smuj - 



(43) 



T"hf,2(w) and Thf,io( w ) are obtained from (25c) and (33). Ta(u;) and Tio(w) have been obtained by numerical 
integration of Q and ([5]). Eq. Q can be reduced to a two-dimensional integral, and the resulting four-dimensional 
numerical integration in Eq. ([5]) was performed carefully to ensure accuracy of O(10 -3 ) for each value of ui. 



E. Dot and Posmom intracules 



-C>HF,2(fc) and -Dhf.ioCO [Eqs. (33) and (25d)], as well as A H f,2(&) and A H F,io(fc) [Eqs. (33) and (25c)] have been 
obtained numerically. Table III gathers the non-interacting and exact (k = 2 and 10) Dot and Posmom intracules in 
Fourier and real space. The similarity between the Dot and Posmom expressions is striking. 



F. Holes 

The correlation hole was originally defined [2] as the difference between the exact and HF Position intracule 

AP(u) = P(u) - Phf(u), (44) 

but this can be extended to any intracule / 

AI = I - J HF . (45) 
One can also define the HF hole as the difference between the HF and non-interacting intracules 



HF 



h 



HF 



(46) 



Fig. [2]shows all of the intracules for n — 2 and Fig. [3]show s the holes created as the Coulomb interaction is introduced. 
One can see from P{u) and M(v) in Figs. 2(a) or| 2(b) that the electrons are found at larger separations and move 



with lower relative momenta in the HF approximation than in the non-interacting case. However, the non-interacting 
and HF intracules, T(cj), D(x) and X(x), are almost identical. The fact that T(cj), D(x) and X(x) are all invariant 
under a uniform scaling leads us to conclude that the introduction of the Coulomb operator at the mean-field level 
leads to an almost exact dilation of the system. 

Fig. [3] reveals that the HF holes and correlation holes of P{u) and M(v) are surprisingly similar in size and shape. 
It also shows that the introduction of correlation decreases the probabilities of co « 7r/2 and x w 0, indicating that 
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FIG. 2. Intracules for a parabolic quantum dot with k = 2: non-interacting( — ), HF ( ) and exact (•■•)• 



the correlated electrons spend less time circularly orbiting their centre of mass. This conclusion is supported by both 
the non-rigorous T(cj) and D(x) intracules and the rigorous X{x) intracule. 



However, the differences between D(x) and X{x) are significant. D(x) comes from the Wigner distribution and, as 
and (16) show, it is the 0(h) approximation to X(x). Its Fourier transform D(k) decays as fc -3 for large 
ihis creates a discontinuity in the second derivative D" (x) at x — [30l 135] . In contrast, X(x) is smooth at 
x = 0. One of the consequences of this misbehaviour at x = is that, for the k = 2 quantum dot, the Dot intracule's 
prediction D(0) — 0.219 overestimates the exact value X(0) = 0.168 by 30%. 



Eqs ([13 

k and t 
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FIG. 3. Intracule holes for the parabolic quantum dots with k = 2 and k = 10: HF hole ( — ) and correlation hole ( ). 
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TABLE IV. HF, radial and exact energies of various helium-like ions. 



Atom 


H 


He 


T -4- 

Li + 


B + 


AT 84- 

Ne s+ 


Z 


1 


2 


3 


5 


10 


HF 


-0.487 93 


-2.86167 


-7.236 41 


-21.986 2 


-93.861 1 


Radial 


-0.5145 


-2.8790 


-7.252 5 


-22.0015 


-93.875 9 


Exact 


-0.52775 


-2.903 72 


-7.279 91 


-22.0310 


-93.906 8 


E c 


-0.03982 


-0.04205 


-0.04350 


-0.0448 


-0.0457 


%E™ d 


66.7 


41.3 


37.0 


34.2 


32.3 



IV. HELIUM-LIKE IONS 



We now turn our attention to the helium-like ions. The Hamiltonian is obtained by substituting the harmonic 



potential V(r) in (27) with the Coulombic potential 



V(r) 



Z 
r 



(47) 



where Z is the nuclear charge. As before, the HF wave function and energy can be found by expanding the HF orbital 
in a Gaussian basis, optimizing both the coefficients and exponents. We consider five values of Z, corresponding to 
the H~, He, Li + , B 3+ and Ne 8+ ions and, in this Section, we focus on their Position and Posmom intracules. 



Ground state 



The HF orbital of the 1 S* ground state was approximated by a Gaussian expansion (33 1 with Nq = 11. The exact 
wave function was approximated by the 64-term Hylleraas-type expansion [50] 



nlm 

We also considered the 64-term radially-correlated wave function [51] [52] 

7 

* rad (ri,r 2 ) = c m(ri + r 2 ) n (ri - r 2 ) 2l e - a ^ +r ^ . 



(48) 



(49) 



Table IV gathers the HF, exact and radially correlated energies obtained from (33), (48) and (49), respectively. 
Only the correct figures are reported [52H54] . as well as the percentage of radial correlation (%££ a ^"r As above, the 



correlation hole is defined as the difference between exact and HF intracules [Eq. (45 )]. We also define the radial and 
angular holes 



A7 r 
AL 



ad 



^rad — ^HF , 



I, 



ad • 



(50) 
(51) 



Both P{u) and X(x) were obtained numerically and are shown in Fig. [4] Increasing the nuclear charge barely 
affects Xi$f{x) because it produces an almost uniform contraction of the system. However, the effect on the exact 
X{x) is much larger and the values of X(0) are 0.1642, 0.1917, 0.1968, 0.2003 and 0.2027 for Z = 1, 2, 3, 5 and 10, 
respectively. These values reveal that the electrons in H~ spend less time mutually orbiting than those in Ne 8+ . This 
is consistent with the conventional view that H _ is a more strongly correlated system than Ne 8+ . 

Fig. [5] shows the Posmom and Position holes of the ions. The depth of AI(x) decreases as Z increases but the 
depth of AP(u) is almost constant as it is squeezed toward the origin. AP(u) exhibits a secondary hole, discussed in 
detail by Pearson et al. [18], but this subtle correlation effect is not visit ble in AX(x). 

Radial correlation provides the majority (67%) of the total correlation energy in H~ but this decreases to 41% 
in He, and to 32% in Ne 8+ , as angular correlation effects becomes dominant. This shift is visible in AF rac i(«) and 
APa,n g {u), as shown in Fig. [m but AX Ta d(x) is always larger than AX ang (x) and becomes almost identical for Ne 8+ . 



If we compare the insets in Fig. |5(b)| and Fig. 6(c)[ we see that the radial secondary hole is several times deeper 
than the total secondary hole. This implies that the radial secondary hole, which has been noted previously by Katriel 
et al. [55], is almost entirely cancelled by an angular secondary hole. 
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FIG. 4. Posmom intracules of helium-like ions: Z = 1 ( ), Z = 2 ( ), Z = 3 (•••), Z = 5 ( ) and Z = 10 (- 
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FIG. 6. Radial and angular holes in He-like ions: Z = 1 (- 



-), Z = 2 (- • -), Z = 3 (• ■ • ), Z = 5 (- - -) and Z = 10 (— ). 
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TABLE V. Origin Posmom intracule, X(0), for various singlet and triplet excited states of the helium atom (i,j = x, y, z) 



Configuration Singlet Triplet 

Is 2 0.2060 

ls2s 0.0866 0.09520 

ls2p 0.0928 0.09525 

2s 2 0.1647 

2s2p 0.1455 0.1153 

2p} 0.1267 

2 Pi 2 Pj 0.1567 0.1538 



X(x) 




FIG. 7. Posmom intracule for states of the He atom: Is 2 X S (— ), ls2s 1 S (---), ls2s 3 5" (• ■ • ), 2s 2 1 S (- ■ -), 2p 2 1 S (- • • -). 



B. Excited States 

We have calculated the HF Posmom intracule for several excited states of the He atom using a Gaussian basis of 
36 s-type functions with exponents 2~ 15 , 2~ 14 , • • • , 2 20 and 31 p-type functions with exponents 2~ 10 , 2~ 14 , ■ • • , 2 20 . 
The maximum overlap method (MOM) has been employed for finding excited-state solutions to the HF self-consistent 
field equations [56] . 

The intracules in the ground state (Is 2 1 S I ) and four excited states (ls2s 1 S, ls2s 3 S, 2s 2 X S and 2p 2 1 S') are shown 
in Fig. [7] Table |v| lists the values of X(0) for these and other excited states. 

When the electron pair occupies a more diffuse orbital, X{x) becomes broader and X(0) drops from 0.206 in the Is 2 
state to 0.165 in the 2s 2 state, and to 0.127 for the 2p 2 state. The decrease is even more marked when the electrons 
occupy orbitals in different shells, such as in the ls2s X S state where -^(0) = 0.087. However, if the two orbitals have 
the same principal quantum number, such as in the 2s2p 1 P state, the decrease is smaller. 

The Dot intracule has been calculated for the first excited state ls2s 3 S [3D] and shows a small dip in D(x) around 
x = 0, which we have previously attributed to the Fermi hole. We now believe that that explanation was incorrect 
and that the dip is a failure of D(x) to capture the behaviour of X(x). 



C. D-dimensional helium atom 



Following the pioneering work of Loeser and Herschbach on the effect of dimensionality on the HF [571 HE] an d 
exact energies |59j of He-like ions, several other two-electron systems have recently been studied in T> dimensions 

H7HBDH53I. 



The generalization of the Posmom intracule for a 2?-dimcnsional space is straightforward. Eqs (10), (11) and (13) 
are unchanged and Eq. (25e) becomes 



JV/2- 



(52) 
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FIG. 8. HF Posmom intracules of the ©-dimensional He atom: V = 2 (— ), V = 3 ( ), 5 = 4 (•••), 5 = 5 (-• -). 



We used a large, even-tempered Gaussian basis, optimizing the coefficients Cj to minimize the HF energy [581 E2 

V 2 



where 



V(r) 



' 'u !•' ('•)'//• -:■ / / '0HF( r i) ^ ^HF( r 2)rfridr 2) 



(53) 



1 



V 2 : 

dr 



3-T> 1 £> min(ri,r 2 ) 
2 ' 2 ' 2 ' max(ri,r 2 ) : 



max(ri, r 2 ) 
d 2 V-l d 
dr 2 r dr ' 



r(x>/2) 



(54) 

(55) 
(56) 



and F is the Gauss hypergeometric function |48j . Our energies for T> = 2, 3, 4 and 5 agree within a microhartree with 
the benchmark values of Herschbach and co-workers [55] • 

Figure [8] shows how X(x) changes with V. The observation that the intracule broadens as T> increases is consistent 
with the conclusion of Herrick and Stillinger [65] that the electrons in X>-helium can avoid each other more easily 
when T> is large. Furthermore, they have shown that the binding energy of the ground state in T> = 5 corresponds 
exactly to the binding energy of the 2pi2pj 3 P state in D = 3. This feature is due to interdimensional degeneracies, 
first noticed by van Vleck [66], and observed for various systems [6T1 l65l 167^70] . We observe likewise that X(x) for 
the Is 2 state in T> = 5 is identical to X(x) for the 2pi2pj 3 P state in T> = 3. 



V. CONCLUSION 



We have introduced a new two-particle density distribution, the Posmom intracule, which condenses information 
about both the relative position and relative momentum of the particles. We have shown how to construct this 
distribution from the many-particle wave function and we have shown that the Dot intracule D{x) in a first-order 
approximation of the Posmom intracule X(x). We have applied our new formalism to two-electron quantum dots and 
the the He-like ions. A comparison between various intracules (Position, Momentum, Angle, Dot and Posmom) has 
been carried out, showing the interrelated information conveyed by these two-particle probability distributions. 
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